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ABSTRACT: This paper presents a summary of the current status of determinations of the 
strong coupling constant a s . A detailed description of the definition, scale dependence and 
inherent theoretical ambiguities is given. The various physical processes that can be used to 
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experimental. 

CONTENTS 



QCD AND ITS COUPLING 



a, FROM Z DECAYS AND e+e" TOTAL RATES 



°This work was supported in part by the Director, Office of Science, Office of High Energy 
and Nuclear Physics, Division of High Energy Physics of the U.S. Department of Energy under 
Contracts DE-AC03-76SF00098. 



DETERMINATION OF a s FROM DEEP INELASTIC SCATTERING 



10 



FROM FRAGMENTATION FUNCTIONS 15 



a-, FROM EVENT SHAPES AND JET COUNTING 15 



y s FROM t DECAY 19 



a, FROM LATTICE GAUGE THEORY COMPUTATIONS 22 



3i s FROM HEAVY QUARK SYSTEMS 23 



FROM HADRON-HADRON SCATTERING 26 



CONCLUSION 27 



2 



3 



1 QCD AND ITS COUPLING 

Quantum chromo dynamics (QCD) is a gauge field theory that describes the 
strong interactions of quarks and gluons Q. All experimental results to date 
are consistent with QCD predictions to within the experimental and theoretical 
errors. In this review, we discuss the current status of the extraction of the strong 
interaction coupling constant a s from the experimental data. 

The QCD Lagrangian describing the interactions of quarks and gluons is 

L = --F^F^ + W ~ m k ) (1) 

k 

where 

F%, = d^Al - d v A° + gf abc A b ^A c u (2) 
is the gluon field strength tensor, 

D^ = d^- igA«T a (3) 

is the gauge covariant derivative, and T a are the 577(3) representation matrices 
normalized so that tiT a T b = 5 ab /2, and the sum on k is over the six different 
flavors (u,d,s,c,b,t) of quarks. At the classical level, the QCD Lagrangian de- 
pends on the six quark masses m^, and the strong interaction coupling constant 
g, or equivalently, the strong fine-structure constant a s = g 2 /A7r. The quantum 
theory contains an additional parameter, the #-angle, that violates CP. The ex- 
perimental limit on this parameter is 6 < 1CP 9 0], so we will set it to zero for 
the purposes of this article. 

One can evaluate QCD scattering amplitudes in powers of a s using a Feynman 
diagram expansion. As is typical in a quantum field theory, loop graphs are 
divergent and need to be treated using a renormalization scheme. The most 
commonly used scheme is modified minimal subtraction (MS) |J, and we will 
use this scheme throughout. An important consequence of renormalization is 
that the parameters a s and m/% of the QCD Lagrangian depend in a calculable 
manner on the MS subtraction-scale \x. The dependence of a s is described by 
the /3-function, 

/iV = 0(«.(M)). (4) 



In perturbation theory, 



0(a s ) = -f3 O7 ^-p l7 ^-p 2 — ^3-..., (5) 



2tt (2vrr (2tt) 

where (for rif flavors of quarks 



O = 11- ~n f , (6) 
19 

Pi = 51- —n f , (7) 



and the next two terms are also known 
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If a s is small, the renormalization group equation Eq. (Q) can be integrated 
using only the (3q term to give 



a s (m) a s (m) 2tt /j, 2 

Since /3o > for nj < 16.5, a s (/i) — > as \x — > oo. The vanishing of the QCD 
coupling for large values of /i is referred to as asymptotic freedom. One important 
consequence of asymptotic freedom, is that QCD processes at high energies can 
be reliably computed in a perturbation expansion in a s . 

A measurable quantity, such as the total cross-section for e + e~ — * hadrons at 
high energies can be computed as a function of the QCD coupling constant a s (/j,) 
and the center of mass energy Ecm, 

a (e + e~ — > hadrons) = — ^ — / (a s (/i) , In — \ . (9) 

E CM V / 

where / is a dimensionless function of its arguments. The form of the cross- 
section given in Eq. (||) follows from dimensional analysis: a has dimensions 
of energy -2 , and /i has dimensions of energy. In the MS sscheme, any depen- 
dence on n is logarithmic, so / can only depend on InE'cM/V- The cross section 
a (e + e~ — > hadrons) is a measurable quantity and cannot depend on the sub- 
traction scale /i, so the [J, dependence on the right hand side of Eq. (|9|) must 
cancel, 

M^/ (a. (M) ,^f) = 0, (10) 

and any value of \i can be used on the right hand side of Eq. (||). In practice, 
one can only compute the right hand side of Eq. (^) at some finite order in 
perturbation theory, and the approximate value of / can depend on n at higher 
order in perturbation theory. Typically, one finds that the perturbation expansion 
has terms of the form 



a s (n) In 



(11) 



which are referred to as "leading logarithms." Even if a s (/i) is small, the per- 
turbation expansion can break down if IhEcm/^ is large. For this reason, it is 
conventional to choose the subtraction scale of order the center of mass energy 
Ecm- The exact choice of scale (for example, whether \x = 2£"cm or -Ecm or 
Ecm/2) is arbitrary, and differences in choice of scale are formally of higher or- 
der in a s . Many methods have been proposed to determine the optimum scale to 
use for a given calculation The only way to determine the "best" scale at a 
given order is to compute the cross-section at next order. [Of course, in this case, 
one might as well use the more accurate formula to determine the cross-section.] 
The scale dependence of a given quantity can also be used to estimate the size 
of neglected higher order corrections. Scale dependence is a dominant source of 
error in many of the quantities that will be used to determine a s . 

Perturbation theory is valid if one chooses /i to be of order Eqm, so that the 
expansion parameter is a s (i?cM), with no large logarithms. This shows that at 
high-energies, the coupling constant is small because of asymptotic freedom, and 
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QCD cross-sections are approximately those of free quarks and gluons. At low, 
energies, non-perturbative effects become important. 

The value of a s is determined by computing a quantity in terms of a s , and 
comparing with its measured value. One might think that it is better to use high- 
energy processes to determine a s , since perturbation theory is more reliable. This 
is not necessarily the case. High energy processes can be computed more reliably 
precisely because they do not depend very much on a s . This means that errors 
in the experimental measurement or theoretical calculation get amplified when 
they are converted to an error on a s . We will see in this article that low-energy 
extractions of a s have comparable errors to those at high energy. 

In addition to scale dependence, the coupling constant a s is also subtraction 
scheme dependent. The scheme dependence of the coupling constant is compen- 
sated for by the scheme dependence in the functional form for a measurable quan- 
tity, so that the value of an observable is scheme independent. The MS scheme 
will be used in this article, but there is still some residual scheme dependence we 
need to consider. In the MS scheme, heavy quarks do not decouple in loop graphs 
at low-energy. For example, the /3-function coefficient (3q = 11 — 2nj/3, where rif 
is the number of quark flavors. This expression is true for all energies, irrespective 
of the mass of the quark. One might expect that at energies much smaller than 
the mass rriQ of a heavy quark, the quark does not contribute to the /3-function. 
This cannot happen in the MS scheme, since MS is a mass-independent sub- 
traction scheme, which results in a mass-independent /3-function. What happens 
is that there are large logarithms of the form In rriQ/fj, that compensate for the 
"incorrect" /3-function at low-energies. In practice, one deals with this problem 
using an effective field theory. At energies smaller than rriQ, one switches from 
a QCD Lagrangian with rif flavors to a QCD Lagrangian with rif — 1 flavors by 
integrating out the heavy quark flavor. The effect of the heavy quark is taken into 
account by higher dimension operators in the QCD Lagrangian, and by shifts in 
the parameters a s and m^. Thus it is necessary to specify the effective theory 
when quoting the value of a s . We will use the notation aT 1 ^ to denote the value 
of a s in the rif flavor theory. One can compute the relation between as and 
a s ni ^ at the scale /i = toq of the heavy quark. This relation is known to three 



When quoting the value of a s (/j,), it is also necessary to specify the number of 
flavors in the effective theory. In most of the a s determinations we consider, the 

(5) 

appropriate effective theory to use is one with rif = 5, and a s will refer to as 
unless otherwise specified. 

Classical QCD is a scale invariant theory, but this scale invariance is broken 
at the quantum level. The quantum theory has a dimensionful parameter A that 



loops H 




(12) 
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characterizes the scale of the strong interactions. The A parameter is determined 
in terms of a s (fx). The solution of the renormalization group equation Eq. (Q) 
including the first three terms in the /^-function is 



47T 



A)ln(// 2 /A 2 ) 



2/3 1 ln[lnQu 2 /A 2 )] 
(31 ln(^ 2 /A 2 ) 



In 



/3 4 ln 2 ( M 2 /A 2 ) \V L \"' J\ 2) Pi 



In fi 2 /A 



+ 



(13) 



This equation can be used to determine A if a s is known at some scale /x. The 
value of A depends on the number of terms retained in Eq. (|l^). The expansion 
parameter in Eq. (O) is 

ln[ln(^ 2 /A 2 )] 
In Gu 2 /A 2 ) ' 1 j 

which is small as long as \i 3> A. In QCD, A is of order 200 MeV. The last term 

in Eq. ( |l3| ) is often dropped in the definition of A. For a fixed value of a s (Mz), 

the shift in A is approximately 15 MeV if the last term in Eq. ( |l3| ) is dropped. 

The QCD /3-function depends on rif, and so changes across quark thresholds. 

This in turn implies that A changes across quark thresholds, so that A^ n f) is 

the value of A with nf dynamical quark flavors. The matching conditions for 
A(n/) A (n 

f ' can be computed using the matching condition Eq. ( |12| ) for a s . 
The differences between A^ , A^\ and A^ 5 ) are numerically very significant. 

In addition to the perturbative effects discussed so far, non-perturbative ef- 
fects play an important role in strong interaction processes. The size of non- 
perturbative effects is governed by the ratio of the strong interaction scale A 
to the typical energy -Ecm of a given process. In many cases, non-perturbative 
effects are estimated using a model analysis, or by a phenomenological fit to 
the experimental data. A few processes can be analyzed rigorously using the 
operator product expansion (OPE) @; in this case, one can calculate the size 
of non-perturbative effects in terms of matrix elements of gauge invariant local 
operators. Two classic examples of this type are deep inelastic scattering, and 
the total cross-section for e + e~ — ► hadrons. 

It is convenient to calculate R(s), the ratio of the total cross-sections for 
e + e~ — > hadrons and e + e~ — > [i + fi~ at center of mass-energy Eqm = Vs- One 
can show that 

R(s) = fo M*)) + h Ms)) {F ^ n + • • • , (15) 

where the first non-perturbative correction depends on the vacuum expectation 
value of the square of the gluon field-strength tensor. By dimensional analysis, 
this quantity is of order A 4 , so the non-perturbative corrections are of order 
A j s 2 . The existence of an OPE provides some crucial information on the size 
of non-perturbative corrections. For R(s), we know that the corrections vanish 
at least as fast as A 4 /s 2 for large values of s, because F^F^ is the lowest 
dimension operator that can contribute in the OPE. Similarly, it is known in 
deep inelastic scattering that the first non-perturbative corrections arise from 
twist-four operators, and are of order A 2 /Q 2 , where Q is the momentum transfer. 
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One can estimate the size of non-perturbative corrections for processes with an 
OPE, by estimating the value of operator matrix elements. In some cases, one 
is fortunate enough that the relevant matrix element can actually be determined 
from some other measurement, or computed from first principles. The size of 
non-perturbative corrections is much less certain if the process does not have an 
OPE. Non-perturbative effects can fall off like a fractional power of A/s, or could 
have some more complicated dependence on s. Typically, one uses some model 
estimate of the non-perturbative corrections. 

Perturbative and non-perturbative corrections to scattering cross-sections are 
interrelated, because the QCD perturbation series is an asymptotic expansion, 
rather than a convergent expansion. A dimensionless quantity f(a s ) has an 
expansion of the form 

f(a s ) = c + cia s + c 2 a 2 s + . . . . (16) 

Typically, the coefficients c n grow as n!, so that the series has zero radius of con- 
vergence. The large-order behavior of the perturbation series can be computed 
in certain limiting cases || & [l0|, 11, 12, If one studies QCD in the limit 



of a large number of flavors, rif — > oo, with a s rif fixed, one can sum all terms 
of the form (ct s nj) n . This is sometimes referred to as the "bubble chain" ap- 
proximation, because the graphs one sums are of the type show in Fig. [l]. QCD 
is not asymptotically free as rif — ► oo. Nevertheless, one can try and apply the 
bubble chain results to QCD. The bubble chain graphs contribute to the QCD 
/3-function. In the large n/ limit, the coefficient (3$ = 11 — 2/3nj —* —2/3nj. One 
therefore computes the bubble chain sum, makes the replacement nj — > — 3(3q/2, 
and uses the resultant expression for QCD with (3q > 0. This seemingly unjus- 
tified procedure has provided some useful insights into the nature of the QCD 
perturbation series. A detailed discussion of this method is beyond the scope of 
the present article. It is typically found that the coefficients c n in the pertur- 
bation expansion have a factorial divergence in the bubble chain approximation. 
One can try and sum a series of this type using a Borel transformation. One 
defines the Borel transform of / by 

f B (t) = c 5(t) +Cl +c 2 t + ... + + . . . . (17) 

n! 

Then the original function can be obtained by the inverse Borel transform, 

/•OO 

f(a s )= / dtfs^e'^. (18) 
Jo 

Suppose the coefficients of f(a s ) have the form 

c n+ i = a n n\, n > 0. (19) 

Then 

OO i 

fs(t) = coS(t) + Yl ^ = coW) + 737^ (20) 

n=0 

and the inverse Borel transform gives 

f(a s ) = co+ dte-^- -. (21) 

Jo 1 — at 



8 



The behavior of the integral is governed by the singularities in the complex t 
plane, which are referred to as renormalons. If a < 0, the integral is well-defined, 
since the singularity at t = 1/a is not along the path of integration. If a > 0, 
the singularity is along the contour of integration. One can regulate the integral 
by deforming the contour around the singularity. The integral depends on the 
precise prescription used. The prescription dependence is related to the pole at 
t = 1/a, which gives a contribution to the integral of the form 





" 1 " 




exp 




= exp 




aa s _ 





2na 



n) (22) 



The value of \i is chosen to be the typical energy in the process, such as the 
momentum transfer Q for deep inelastic scattering. The perturbative series then 
has the same structure as a typical non-perturbative correction, a power law 
correction of the form 



'A 2 V° 



Q 



2 



(23) 



where no = /?o/(4vra) is the renormalon singularity in the variable u = 0ot/(4ir). 
[It is conventional to refer to the location of the renormalon singularity in u rather 
than in t.j Contributions of the form Eq. ( p3[ ) are called renormalon ambiguities, 
since their value depends on the way in which one performs the inverse Borel 
transform. 

Renormalon ambiguities have the same structure as non-perturbative correc- 



tions. It has been suggested in the literature [13] that renormalons can be 
used as a guide to the size of non-perturbative corrections. There is one non- 
trivial check to this idea. A renormalon singularity at Uq corresponds to a non- 



perturbative ambiguity of the form Eq. (23). In processes that have an OPE, all 



non-perturbative effects should be given by the matrix elements of gauge invariant 
local operators. To every renormalon ambiguity in the perturbation expansion, 
there should be a corresponding ambiguity in the operator matrix element, such 
that the sum is well-defined juj. This can only happen provided that there is a 
gauge invariant local operator corresponding to every renormalon singularity. For 
example, the first gauge invariant operator corrections to deep inelastic scattering 
are of the form A 2 /Q 2 , and it is known that the first renormalon ambiguity is 
at u = 1. Similarly, for R, the first non-perturbative corrections are of the form 
A 4 /Q 4 , and the first renormalon ambiguity is at u = 2. The matching between 
renormalon singularities and the OPE occurs in all examples that have been com- 
puted so far. For this reason, renormalon singularities have also been taken as an 
indication of the size of non-perturbative effects in processes without an OPE. 
Non-perturbative effects are expected to fall off faster if the renormalons are at 
larger values of u. 

Infrared sensitivity is also used to estimate the size of non-perturbative correc- 
tions to a measurable quantity [O]. One computes the quantity in the presence 
of an infrared cutoff momentum A. For example, one can imagine working in 
a box of size 1/A, or using a gluon mass of order A. Cross-sections can have 
infrared divergences of the form In A. If one computes measurable cross-sections 
for color singlet states to scatter into color singlet states, one finds that the In A 
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terms cancel, and the cross-section is infrared finite, a result known as the KLN 



14] theorem. It is important to include finite detector resolution to get a finite 



cross-section, as for QED. For example, the Bhabha scattering cross-section for 
e + e~ — ► e + e~ has an infrared divergence at one-loop order. However, it is im- 
possible to distinguish e + e~ — > e + e~ from e + e~ — > e + e~j if the photon energy 
E-y is smaller than the detector resolution S. The measurable quantity is the 
sum of the e + e~ — > e + e~ cross-section and the e + e~ — ► e + e~ r y cross-section for 
Ely < 5, which is free of infrared singularities. While the In A term must cancel, 
terms of order A, A 2 In A, etc. which vanish as A — > need not cancel. The first 
non- vanishing term is an indication of the infrared-sensitivity of a given quantity. 
In QCD, one can imagine that the scale A represents the confinement scale A. 
A process that is infrared sensitive at order A ra would then be expected to have 
non-perturbative corrections of order A n /Q n , where Q is the typical momentum 
transfer. In a few cases, one can analyze the problem using renormalon methods, 
and by using the criterion of infrared sensitivity. It is found in these cases that 
both methods give the same estimate for the size of non-perturbative corrections. 



2 a, FROM Z DECAYS AND 



e ' e 



TOTAL RATES 



The total cross section for e + e~ — > hadrons is obtained (at low values of y/s) by 
multiplying the muon-pair cross section by the factor R. At lowest order in QCD 
perturbation theory R = R° = 3S g eg where e q is the electric charge of the quark 
of flavor q. The higher-order QCD corrections to this are known, and the results 
can be expressed in terms of the factor: 



R = R<f>) 



IT \ 7T 



a. 



7T 



+ 



(24) 



where C 2 = 1.411 and C 3 = -12.8 p]. 

This result is only correct in the zero-quark-mass limit. The 0(a s ) corrections 
are also known for massive quarks [|l7]] . The principal advantages of determining 
a s from R in e + e~ annihilation are that the measurement is inclusive, that there is 
no dependence on the details of the hadronic final state and that non-perturbative 
corrections are suppressed by 1/s 2 . 

A measurement by CLEO @ at y/s = 10.52 GeV yields q s (10.52 GeV) = 
0.20 ± 0.01 ± 0.06 which corresponds to a s (M z ) = 0.13 ± 0.005 ± 0.03. A com- 



parison of the theoretical prediction of Eqn. 24 (corrected for the 6-quark mass), 
with all the available data at values of y/s between 20 and 65 GeV, gives [fl9| ] 
a,, (35 GeV) = 0.146 ±0.030. It should be noted that the size of the order term 
is of order 40% of that of the order a 2 and 3% of the order a s . If the order 
term is not included, the extracted value decreases to a s (35 GeV) = 0.142±0.03, 
a difference smaller than the experimental error. 

Measurements of the ratio of the hadronic to leptonic width of the Z at 
LEP and SLC, T^/T^ probe the same quantity as R. Using the average of 
lyi^ = 20.783 ± 0.029 gives a s {M z ) = 0.123 ± 0.004 @. The prediction 
depends upon the couplings of the quarks and leptons to the Z. The preci- 
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sion is such that higher order electroweak corrections to these couplings must 
be included. There are theoretical errors arising from the values of top-quark 
and Higgs masses which enter in these radiative corrections. Hence, while this 
method has small theoretical uncertainties from QCD itself, it relies sensitively on 
the electroweak couplings of the Z to quarks J2l] and on the ability of the Stan- 
dard Model of electroweak interactions to predict these correctly. The presence 
of new physics which changes these couplings via electroweak radiative correc- 
tions would invalidate the extracted value of a s (Mz)- Since the Standard Model 
fits the measured Z properties well, this concern is ameliorated and more precise 
value of ot s can be obtained by using a global fit to the many precisely measured 



properties of the Z boson and the measured W and top masses. This gives |22| 

a s (M z ) = 0.1192 ±0.0028 

This error is larger than the shift in the value of a s (Mz) (~ 0.002) that would 
result if the order a s (Mz) 3 term were omitted and hence one can conclude that 
it is very unlikely that the uncertainty due to the unknown a s (Mz) 4: terms will 
dominate over the experimental uncertainty. 



3 DETERMINATION OF a s FROM DEEP INELASTIC SCAT- 
TERING 

The original and still one of the most powerful quantitative tests of perturbative 
QCD is the breaking of Bjorken scaling in deep-inelastic lepton-hadron scattering. 
Consider the case of electron-proton scattering (ep — > eX), where the cross- 
section can be written as 



da iira 2 



em' 



dxdy Q 4 



' + ! ' y) -2xFx(x,Q 2 ) + (1 - y)(F 2 (x,Q 2 ) - 2xF 1 (x,Q 2 )) 



2 

(25) 

The variables are defined as follows (see Figure |2j): q is the momentum of the 
exchanged photon, P is the momentum of the target proton, k is that of the 
incoming electron, and 



m p 

Q 2 

2m p u 

q-p 

k-p 

,2 



2p ■ k + m v 

(26) 



For charged current scattering, which proceeds via the exchange of a virtual 
W boson between the lepton and target nucleus, there is an additional parity 
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violating structure function F< 



3 



A„vN /~i2 7if4 _ 



+(i-i,-,yM 2 /Q 2 )Ff(^ 2 ) ( 27 ) 

-^((i-^-i^^g 2 )) 

For VN scattering the sign of the last (xF^) term is reversed. 

In the leading-logarithm approximation, the measured structure functions Fi(x, Q 2 ) 
are related to the quark distribution functions qi(x,Q 2 ) according to the naive 
parton model, for example 

F 2 (x,Q 2 ) = £>?<Zi(x,Q 2 ) (28) 

Here qi{x, Q 2 ) is the probability for a parton of type i to carry a fraction x of the 
nucleon's momentum. The Q 2 dependence of the parton distribution functions 
is predicted by perturbative QCD, hence a measurement of the Q 2 dependence 
( "scaling violation" ) can by used to measure a s . In describing the way in which 
scaling is broken in QCD, it is convenient to define nonsinglet and singlet quark 
distributions: 

F NS = qi - qj = 

i 

The nonsinglet structure functions have nonzero values of flavor quantum num- 
bers such as isospin or baryon number. The variation with Q 2 of these is described 



by the so-called DGLAP equations |2J J2J : 



2 dF NS _ a s (\Q\) pqq rpNS 

~dQ 2 ~ ~ 2tt r * P 



Q 2 9 



as(\Q\) ( pqq 2n f pq3 
2tt \ 



pai pgg 



W 2 \G ) 

where * denotes a convolution integral: 

f*9= [ — f(y) 9 (- 
Jx y \y 



Co) 



The leading-order Altarelli-Parisi [24] splitting functions are 



pqq 
pqg 
pgq 
pgg 



l+x 2 



+ 25(1 - x) 



{l-x) 

\ [J + (1 - x) 2 ] 



l+(l-x) 2 

x 

-%6(l-x) 
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Here the gluon distribution G(x, Q 2 ) has been introduced and 1/(1 — x)+ means 



o (1 — a;)+ Jo 

Measurement of the structure functions over a large range of x and Q 2 allows both 
a s and the parton distributions to be determined. Notice that a s and the gluon 
distribution can only be obtained by measuring the Q 2 dependence. The precision 
of contemporary experimental data demands that higher-order corrections also 



be included [25] . The above results are for massless quarks. Algorithms exist for 
the inclusion of nonzero quark masses p6[ . These are particularly important for 
neutrino scattering near the charm threshold. At low Q 2 values, there are also 
important "higher-twist" (HT) contributions of the form: 

FlM 2 )= F r ) (^Q 2 ) + FfT Q ( 2 x ' Q2) + --- 

Leading twist (LT) terms are those whose behavior can be predicted using the 
parton model, and are related to the parton distribution functions. Higher-twist 
corrections depend on matrix elements of higher dimension operators. These 
corrections are numerically important only for Q 2 <C(few GeV 2 ) except for x 
very close to 1. At very large values of x corrections proportional to log(l — x) 



can become important [27] 



From Eqn. |3[ it is clear that a nonsinglet structure function offers in principle 
the most precise test of the theory, since the Q 2 evolution is independent of 
the unmeasured gluon distribution. The CCFR collaboration fit to the Gross- 
Llewellyn Smith sum rule ^] which is known to order ||29|| |30f| (Estimates of 
the order a* term are available Q]) 

J 1 dx{F^(x,Q 2 ) + F^(x,Q 2 )} = 

3 [1 - + 3.58^ + 19.0(^) 2 )] - AHT 
where the higher-twist contribution AHT is estimated to be (0.09 ± 0.045)/Q 2 in 



and to be somewhat smaller by [33]. The CCFR collaboration |3iJ], com 



bines their data with that from other experiments |35|] and gives a s (\/3 GeV) 
0.28 ±0.035 (expt.)±0.05 (sys)l^| 5 (theory). The error from higher-twist terms 
(assumed to be AHT = 0.05±0.05) dominates the theoretical error. If the higher 



twist result of [33 is used, the central value increases to 0.31 in agreement with 
the fit of ||]. This value extrapolates to a s (M z ) = 0.118 ± 0.011. 

Measurements involving singlet-dominated structure functions, such as F2, re- 
sult in correlated measurements of a s and the gluon structure function. A full 
next to leading order fit combining data from SLAC @, BCDMS §8|, E665 g 



and HERA |4(|has been performed |41j| . These authors extend the analysis to 
next to next to leading order (NNLO). In this case the full theoretical calcula- 
tion is not available as not all the three-loop anomalous dimensions are known; 
their analysis uses moments of structure functions^ and is restricted to those mo 

The moments are denned by M n = / Q x n F{x,Q 2 )dx. 
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ments where the full calculation is available [42], 36]. The NNLO result is 
a s (M z ) = 0.1172 ± 0.0017 (expt.) ± 0.0017 (sys). Here the first error is a com- 
bination of statistical and systematic experimental errors, and the second error 
is due to the uncertainties, quark masses, higher twist and target mass correc- 
tions, and errors from the gluon distribution. If only a next to leading order fit 
is performed then the value decreases to a s (M z ) = 0.116 indicating that the 
theoretical results are stable. No error is included from the choice of \i\ /i = Q 
is assumed. We use a total error of ±0.0045 to take into account an estimate of 
the scale uncertainty. This result is consistent with earlier determinations ||43J, 



|4j, and @. 

The spin-dependent structure functions, measured in polarized lepton nucleon 
scattering, can also be used to determine a s . The spin structure functions G\ and 
(j?2 are defined in terms of the asymmetry in polarized lepton nucleon scattering 

^ ' ^ dxdy dxdy ^ ^ 

where the subscript p (ap) refers to the state where the nucleon spin is parallel 
(anti-parallel) to its direction of motion in the center of mass frame of the lepton- 
nucleon system. In both cases the lepton has its spin aligned along its direction 
of motion. 



''<■'••//> - ^((l-2/ ! / 2 + 2xVM 2 /Q 2 )G 1 (,,Q 2 ) 

MQ 2 ( 30) 

+4x 2 M 2 G 2 (x,Q 2 )/Q 2 ) 

The Q 2 evolution of the spin structure functions Gi(x,Q 2 ) and G2(x,Q 2 ) is 
similar to that of the unpolarized ones and is known at next to leading order 
1 51]. Here the values of Q 2 ~ 2.5 GeV 2 are small particularly for the E143 data 



1 46] and higher-twist corrections are important. A fit |47j using the measured spin 
dependent structure functions measured by themselves and by other experiments 
P§ gives a s (M z ) = 0.121 ± 0.002(expt.) ± 0.006(theory and syst.). Data 



from HERMES JS] are not included in this fit; they are consistent with the 
older data. 



a s can also be determined from the Bjorken sum rule [50]. 



5 Bj = C dx (G\ - Gi) = -U (31) 
Jo o 

At lowest order in QCD a 3 = g A = = 1.2573 ± 0.0028. A fit gives || 
a s (Mz) = 0.118^o;q24; a significant contribution to the error being due to the 
extrapolation into the (unmeasured) small x region. Theoretically, the sum rule 
is preferable as the perturbative QCD result is known to higher order, and these 
terms are important at the low Q 2 involved. It has been shown that the theoret- 
ical errors associated with the choice of scale are considerably reduced by the use 
of Pade approximants |53| which results in a s (1.7 GeV) = 0.328 ± 0.03(expt.) ± 
0.025(theory) corresponding to a s {M z ) = 0.116±^| (expt.)±0.003(theory). No 
error is included from the extrapolation into the region of x that is unmeasured. 
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Should data become available at smaller values of x so that this extrapolation 
could be more tightly constrained, the sum rule method could provide a better 
determination of a s than that from the spin structure functions themselves. 

At very small values of x and Q 2 , both the x and Q 2 dependence of the structure 
functions is predicted by perturbative QCD |54j. Here terms to all orders in 
a s ln(l/x) are summed. The data from HERA on F^ix, Q 2 ) can be fitted to 
this form |^], including the NLO terms which are required to fix the Q 2 scale. 
The data are dominated by 4 GeV 2 < Q 2 < 100 GeV 2 . The fit @] using HI 
data gives a s (M z ) = 0.122 ±0.004 (expt.) ± 0.009 (theory). (The theoretical 
error is taken from |{55f .) The dominant part of the theoretical error is from the 
scale dependence; errors from terms that are suppressed by l/log(l/x) in the 
quark sector are included [58] while those from the gluon sector are not. 

Typically, A is extracted from the deep inelastic scattering data by parameter- 
izing the parton densities in a simple analytic way at some Qq, evolving to higher 
Q 2 using the next-to-leading-order evolution equations, and fitting globally to 
the measured structure functions. Thus, an important by-product of such stud- 
ies is the extraction of parton densities at a fixed-reference value of Qq. These 
can then be evolved in Q 2 and used as input for phenomenological studies in 
hadron-hadron collisions (see below). These densities will have errors associated 
with the that value of a s . A next-to-leading order fit must be used if the process 
being calculated is known to next-to-leading order in QCD perturbation theory. 
In such a case, there is an additional scheme dependence; this scheme depen- 
dence is reflected in the 0{a s ) corrections that appear in the relations between 
the structure functions and the quark distribution functions. There are two com- 
mon schemes: a deep-inelastic scheme where there are no order a s corrections in 
the formula for F2(x,Q 2 ) and the minimal subtraction scheme. It is important 
when these next-to- leading order fits are used in other processes (see below) , that 
the same scheme is used in the calculation of the partonic rates. Most current 
sets of parton distributions are obtained using fits to all relevant data | |59| |. In 
particular, data from purely hadronic initial states are used as they can provide 
important constraints on the gluon distributions. 



3.0.1 PHOTON STRUCTURE FUNCTIONS 

Experiments in e + e~ collisions can be used to study photon-photon interactions 
and to measure the structure function of a photon |6(| , by selecting events of the 
type e + e~ — ► e + e~ +hadrons which proceeds via two photon scattering. If events 
are selected where one of the photons is almost on mass shell and the other has 
a large invariant mass Q, then the latter probes the photon structure function 
at scale Q; the process is analogous to deep inelastic scattering where a highly 
virtual photon is used to probe the proton structure. The Q 2 variation of this 
structure function follows that shown above (see Eq |3|) . 

A review of the data can be found in |6l| . Data have become available from 
LEP [32] and from TRISTAN [^] [64] which extend the range of Q 2 to of order 



300 GeV 2 and x as low as 2 x 10 3 and show Q 2 dependence of the structure 
function that is consistent with QCD expectations. Experiments at HERA can 
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also probe the photon structure function by looking at jet production in 7p col- 
lisions; this is analogous to the jet production in hadron-hadron collisions which 
is sensitive to hadron structure functions. The data p!| are consistent with 
theoretical models l66|]. 



a s FROM FRAGMENTATION FUNCTIONS 



Measurements of the fragmentation function di(z,E), the probability that a 
hadron of type i be produced with energy zE in e + e~ collisions at y/s = 2E, 
can be used to determine a s . As in the case of scaling violations in structure 
functions, QCD predicts only the E dependence in a form similar to the Q 2 de- 
pendence of Eq[3|. Hence, measurements at different energies are needed to extract 
a value of a s . Because the QCD evolution mixes the fragmentation functions for 
each quark flavor with the gluon fragmentation function, it is necessary to deter- 
mine each of these before a s can be extracted. The ALEPH collaboration has 
used data in the energy range yfs = 22 GeV to = 91 GeV. A flavor tag is used 
to discriminate between different quark species, and the longitudinal and trans- 
verse cross sections are used to extract the gluon fragmentation function ||67|] . 
The result obtained is a s (M z ) = 0.126 ± 0.007 (expt.) ± 0.006 (theory) |jf]. 
The theory error is due mainly to the choice of scale at which a s is evaluated. 
The OPAL collaboration |69] has also extracted the separate fragmentation func- 



tions. DELPHI [7C] has performed a similar analysis using data from other 
experiments at center of mass energies between 14 and 91 GeV with the result 
a s (Mz) = 0.124 ± 0.007 ± 0.009 (theory). The larger theoretical error is because 
the value of /i was allowed to vary between 0.5-y/i and 2yfs. These results can be 
combined to give a s (M z ) = 0.125 ± 0.005 ± 0.008 (theory). 



5 a s FROM EVENT SHAPES AND JET COUNTING 

An alternative method of determining a s in e + e~ annihilation involves measuring 
the the topology of the hadronic final states. There are many possible choices 
of inclusive event shape variables: thrust [[n]], energy-energy correlations [fr^ |, 
average jet mass, etc.. These quantities must be infrared safe, which means that 
they are insensitive to the low energy properties of QCD and can therefore be 
reliably calculated in perturbation theory. For example, the thrust distribution 
is defined by 

T = rnax(^2\pi-n\/^2\pl\), (32) 

i i 

where the sum runs over all hadrons in the final state and the unit vector n is 
varied. At lowest order in QCD the process e + e~ — > qq results in a final state 
with back to back quarks i.e. "pencil- like" event with T = 1. Alternatively, the 
event can be divided by a plane normal to the thrust axis and the invariant mass 
of the particles in the two hemispheres is computed, the larger (smaller) of these 
is M h (Mi). At lowest order in QCD M h = M t = 0. 
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The observed final state consists of hadrons rather than the quarks and gluons 
of perturbation theory. The hadronization of the partonic final state has an 
energy scale of order A. The resulting hadrons acquire momentum components 
perpendicular to the original quark direction of order A. This effect induces 
corrections to the shape variables of order A/y/s. A model is needed to describe 
the detailed evolution of a partonic final state into one involving hadrons, so that 
detector corrections can be applied. Furthermore if the QCD matrix elements 
are combined with a parton-fragmentation model, this model can then be used to 
correct the data for a direct comparison with the perturbative QCD calculation. 
The different hadronization models that are used [^] model the dynamics that 
are controlled by non-perturbative QCD effects which we cannot yet calculate. 
The fragmentation parameters of these Monte Carlo simulations are tuned to get 
agreement with the observed data. The differences between these models can be 
used to estimate systematic errors. 

In addition to using a shape variable, one can perform a jet counting experi- 
ment. At order a s the partonic final state qqg appears which can manifest itself 
as a three-jet final state after hadronization. Every higher order produce a higher 
jet multiplicity and measuring quantities that are sensitive to the relative rates 
of two-, three-, and four-jet events can lead to a determination of a s . There 
are theoretical ambiguities in the way that particles are combined to form jets. 
Quarks and gluons are massless, whereas the observed hadrons are not, so that 
the massive jets that result from combining them cannot be compared directly 
to the massless jets of perturbative QCD. 

The jet-counting algorithm, originally introduced by the JADE collaboration 74], 



has been used by many other groups. Here, particles of momenta p, and pj are 
combined into a pseudo-particle of momentum p, + pj if the invariant mass of 
the pair is less than yoy/s. The process is then iterated until no more pairs of 
particles or pseudo-particles remain. The remaining number of pseudo-particles 
is then defined to be the number of jets in the event, and can be compared 
to the perturbative QCD prediction which depends on yo- The Durham algo- 
rithm is slightly different: in computing the mass of a pair of partons, it uses 
M 2 = 2mm(E 2 ,E%)(l — cos 9ij) for partons of energies E{ and Ej separated by 
angle 9y [fT5[ . Different recombination schemes have been tried, for example com- 
bining 3-momenta and then rescaling the energy of the cluster so that it remains 
massless. These varying schemes result in the same data giving slightly different 
values ]7(| \TJ_\ of a s . These differences can be used to estimate a systematic error. 



However, such an error may be conservative as it is not based on a systematic 
approximation. 

The starting point for all these quantities is the multijet cross section. For 
example, at order a s , for the process e + e~ — > qqg: pl| 

1 d 2 a 2a s x\ + x\ 



a dx\dx2 3ir (1 — x{)(l — x-i) 

where Xi = ^-are the center-of-mass energy fractions of the final-state (massless) 
quarks. The order a 2 corrections to this process have been computed, as well as 



the 4-jet final states such as e + e — * qqgg [82]. A distribution in a "three-jet" 
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variable, such as those listed above, is obtained by integrating this differential 
cross section over an appropriate phase space region for a fixed value of the 
variable. Thus < 1 - T >~ a s , < M\ > /s ~ a s and < Mf > /s ~ a 2 s . 

The result of this integration depends explicitly on a s but scale fi at which 
a s {fi) is to be evaluated is not clear. In the case of jet counting, the invariant 
mass of a typical jet (or ^/syo) is probably a more appropriate choice than the 
e + e~ center-of-mass energy. While there is no justification for doing so, if the 
value of /x is allowed to float in the fit to the data, the fit improves and the 
data tend to prefer values of order -y/s/10 GeV for some variables [77| Q]; the 
exact value depends on the variable that is fitted. Typically experiments assign 
a systematic error from the choice of \x by varying it by a factor of 2 around the 
value determined by the fit. The choice of this factor is arbitrary 

Estimates for the non-perturbative corrections to < 1 — T > have been made 
1 84] using an operator product expansion. 

< 1 - T >= A^- + B{^f + (33) 

2-7T 27T y/S 

where A and B known quantities |82| , [i is the renormalization scale and cto is the 
non-perturbative parameter (the matrix element of an appropriate operator) to 
be determined from experiment. Note that the corrections are only suppressed 
by -y/s. This provides an alternative to the use of hadronization models for 
estimating these non-perturbative corrections. The DELPHI collaboration [^] 
uses data below the Z mass from many experiments and Eq. 33 to determine 
a s (Mz) = 0.119 ± 0.006, the error being dominated by the choice of scale. The 
values of a s and the non-perturbative parameter ao are also determined by a 
fit to using the variable < M% > /s. While the extracted values of a s (Mz) are 
consistent with each other, the values of ao are not. The analysis is useful as one 
can directly determine the size of the 1/E corrections; they are approximately 
20% (50%) of the perturbative result at ^/s = 91(11) GeV. Even at *Js = 91 GeV 
the omission of these perturbative terms will cause a shift on the extracted value 
of a s of ~ 0.05 which is much larger than the quoted experimental errors. 

The perturbative QCD formulae can break down in special kinematical con- 
figurations. For example, the first term in Eq. ^ contains a term of the type 
a s ln 2 (l — T). The higher orders in the perturbation expansion contain terms of 
order a™ln m (l — T). For T ~ 1 (the region populated by 2-jet events), the per- 
turbation expansion in a s is unreliable. The terms with n < m can be summed 
to all orders in a s |8(|. If the jet recombination methods are used, higher-order 
terms involve a"ln m (yo), these too can be resummed [jSTj]. The resummed re- 
sults give better agreement with the data at large values of T. Some caution 
should be exercised in using these resummed results because of the possibility of 
overcounting; the showering Monte Carlos that are used for the fragmentation 
corrections also generate some of these leading-log corrections. Different schemes 
for combining the order o? s and the resummations are available [|88| . These differ- 
ent schemes result in shifts in a s (Mz) of order ±0.002. The use of the resummed 
results improves the agreement between the data and the theory. 

Studies on event shapes have been undertaken at lower energies at TRIS- 
TAN, PEP/PETRA, and CLEO. A combined result from various shape pa- 
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rameters by the TOPAZ collaboration gives a s (58 GeV) = 0.125 ± 0.009, us- 
ing the fixed order QCD result, and a s (58 GeV) = 0.132 ± 0.008 (correspond- 
ing to a s (Mz) = 0.123 db 0.007) where the error is dominated by scale and 
fragmentation uncertainties. The CLEO collaboration fits to the order a 2 re- 
sults for the two jet fraction at ^/s = 10.53 GeV, and obtains a s (10.93) = 
0.164 ±0.004 (expt.) ±0.014 (theory) Q. The dominant systematic error arises 
from the choice of scale (//) , and is determined from the range of a s that results 
from fit with \i = 10.53 GeV, and a fit where \i is allowed to vary to get the 
lowest x 2 - The latter results in fi = 1.2 GeV. Since the quoted result corresponds 
to a s (1.2) = 0.35, it is by no means clear that the perturbative QCD expression 
is reliable and the resulting error should, therefore, be treated with caution. A 
fit to many different variables as is done in the LEP/SLC analyses would give 
added confidence to the quoted error. 

Recently studies have been carried out at energies between ~130 GeV |)l]] and 
~189 GeV [||. These can be combined to give a s (130 GeV) = 0.114 ± 0.008 
and a s (189 GeV) = 0.1104 ± 0.005. The dominant errors are theoretical and 
systematic and, as most of these are in common at the different energies, these 
data, those at the Z resonance and lower energy provide very clear confirmation 
of the expected decrease in a s as the energy is increased. 

A combined analysis of the data between 35 and 189 GeV using data from 
OPAL and JADE using a large set of shape variables shows excellent agree- 
ment with a s (M z ) = 0.1187±g;ggf|. A comparison of this result with those at the 
Z resonance from SLD |77|], OPAL f7|, L3 f7|, ALEPH g], and DELPHI f9|, 
indicates that they are all consistent with this value. The experimental errors are 
smaller than the theoretical ones arising from choice of scale \x and modeling of 
non-perturbative effects, which are common to all of the experiments. The SLD 
collaboration |77|] determines the allowed range of [i by allowing any value that is 
consistent with the fit. This leads to a larger error (~ 0.0056) than that obtained 
by DELPHI who vary fx by a factor of 2 around the best fit value and obtain 
±0.0008. We elect to use a more conservative average of a s {Mz) = 0.119±0.005. 

At lowest order in a s , the ep — ► e X scattering process produces a final state of 
(1+1) jets, one from the proton fragment and the other from the quark knocked 
out by the underlying process e + quark — ► e + quark. At next order in a s , a 
gluon can be radiated, and hence a (2+1) jet final state produced. By comparing 
the rates for these (1+1) and (2+1) jet processes, a value of a s can be obtained. 
A NLO QCD calculation is available The basic methodology is similar to 

that used in the jet counting experiments in e + e~ annihilation discussed above. 
Unlike those measurements, the ones in ep scattering are not at a fixed value of 
Q 2 . In addition to the systematic errors associated with the jet definitions, there 
are additional ones since the structure functions enter into the rate calculations. 
Results from HI |97| and ZEUS @ can be combined to give a s (M z ) = 0.118 ± 
0.0015 (stat.) ± 0.009 (syst.). The contributions to the systematic errors from 
experimental effects (mainly the hadronic energy scale of the calorimeter) are 
comparable to the theoretical ones arising from scale choice, structure functions, 
and jet definitions. The theoretical errors are common to the two measurements; 
therefore, we have not reduced the systematic error after forming the average. 
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6 a s FROM r DECAY 

The coupling constant a s can be determined from an analysis of hadronic r 
decays [99, |100] , 101 ]. The quantity that will be used is the ratio 

r(r -> t/ T + hadrons + (7)) 

-"r = =r, = — : , y. , (34) 

r(r -> ^ T ei/ e + (7)) 

where (7) represents possible electromagnetic radiation, or lepton pairs. In the 
absence of radiative corrections, the ratio R T is 

R T = 3 (\V ud \ 2 + \V US \ 2 ) «3, (35) 



where 3 is the number of colors. The experimental value R T = 3.61 ±0.05 is close 
to three, which is experimental evidence for the existence of three colors in QCD. 
The deviation of R T from three is used to extract a s . 
The weak decay Lagrangian for non-leptonic r decay is 

L = -^/f C r(ri [V: d Ur^P L Td^P L U + V^V^ P L TS^P L u] , (36) 

where V us and V ud are the CKM mixing angles. The Lagrangian Eq. ( |3"6| ) is ob- 
tained at the scale fi = Mw by integrating out the W boson to generate a local 
four-Fermion operator in the effective theory below M\y , and C T = 1 at fx = M\y . 
The typical momentum transfer in r decays is of order m T , so it is necessary to 
scale the Lagrangian Eq. (|36|) from fi = My/ to fj, = m T . Electromagnetic interac- 
tions renormalize the Lagrangian. At one-loop, the renormalization from graphs 
shown in Fig. ([|) produce a multiplicative renormalization of the Lagrangian, 
and give 

C T (m T ) = l + ^ln^« 1.009. (37) 

The r decay amplitude r — ► v T X, where X is the final hadronic state, can be 
written as 

A = -i-Z.C T ^)u{j> v )^P L u{p T ) [V; d {X\ d llx P L u |0) + V: S (X\ S^Plu |0)] . 

(38) 

Squaring the amplitude, and computing the decay rate gives 

p _ , Gf \ C r\ 2 \Vud\ 2 v- f d 3 p u r • r> n^ u ie^ 

1-4 ^ ^ J ^f2E v [PuPt+PuPt Pv Pt9 ie a pP v P T 

x{2Ttf5\p T -p v - p x ){0\ u lv P L d \X){X\ d ltl P L u |0), (39) 

where we have retained only the V ud term for simplicity. The sum on X is 
symbolic for the sum over all final states, including phase space factors. The 5 
function can be written as 

& A {Pt -Pu-Px)= J d 4 q5 4 (p T -p„- q)8 A {q - p x ), (40) 

and the sum on X can be written as 

Y,(2*) i S' l (q-Px)(0\u 1 vPLd\X}(X\d^P LU \0) = wM(q). (41) 
x 
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The tensor Wn V is related to another quantity !!„„, defined by 



4^(9) 



d qe 



-iq-x 



T{jl{x)jM |0), 



(42) 



where j 7 ' = dj^Piu. Inserting a complete set of states in the time-ordered 
product, one finds that 

iy« = 2Imn« (43) 
The tensor depends on the only variable, q, and must have the form 



by Lorentz invariance. The tensor is then given by 

W^ d \q) = (-q 2 9 „ + q,q v ) 4"V) + q,q^ d \q") 

where 



n 



(ud) 



2 Imn 



(ud) 



n 



(ud) 



2 Imn 



(ud) 
T ■ 



(44) 



(45) 



(46) 



If the light quark mass difference m^— m u and m s —m u are neglected, the hadronic 



currents d^ u Pi,u and s^ u Plu are conserved. This implies that q^Ii U v 
U L (q 2 ) = 0. Inserting Eq. @) and Eq. (||) into Eq. (|47]) gives 

2 



0, so that 



r 



^Gp \ C T \ 

2 „2> 



(m 



where 



n T (m 2 T + 2q 2 )+n L m 2 T 
2nW 1 HA I 2 q( us ) 



(2^) 3 2^ 



(47) 
(48) 



^T,L = \ Vud\ ^t,L + 1^*1 **t,L 

and we added back the V us contribution. There is no interference term (at lowest 
order in the weak interactions), because the u — > d and u — > s currents lead to 
final states with different flavor quantum numbers, and the strong interactions 
conserve flavor. The hadronic invariant mass distribution can then be written as 

2 



dr 

ds 



Gp \ C T 
8vr 2 



-(m T 



n T (s){m 2 T + 2s) + Q L {s)m< 7 



The ratio of the hadronic to leptonic decay rate of the r is given by [10C] 



R T ' ' ' » K V j 1 I T -^(l- — 



mi 



mi 



9, rl «)[i + ?L)+n L (s 



(49) 



(50) 



The hadronic tensors n^T^) are analytic in the complex s plane, except for 
a branch cut along the positive real axis. The discontinuity across the cut is 
Ql,t(s), and is the cross-section for the currents to create hadrons. Clearly, 
the hadron production rate is sensitive to non-perturbative effects, and can not 
be computed reliably. Far away from the physical cut, there are no infrared 
singularities, and QCD perturbation theory is valid. One can rewrite the integral 
Eq. pi) as 



Rt 



6iri 



ds 

C VPb% 



1 



mt 



n T (s) 1 + 



2s 



nit 



+ n L ( s ) 



(51) 
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where C = C\ is the contour shown in Fig. (||). The difference of II^t above and 
below the cut is SIl,t, so this gives back Eq. (|50[). Since Hl,t have no singularities 
in the complex s plane other than the branch cut, one can deform the contour 
C\ to the contour C2, and use Eq. (pTj) with C = Ci- The advantage of using 
Eq. ( |5l| ) with C = C2 rather than C = C\ is that one needs to know U^t far 
away from the cut for most of the integration contour. The contour approaches 
the cut at s = m 2 ., but at this point, the integrand vanishes as (s — m T ) 2 , so the 
contribution of the region near s = m 2 . to the total integral is suppressed | ]100| | . 

II(s) can be computed in perturbation theory using an operator product ex- 
pansion, which is valid away from the physical cut. The perturbation theory 
result for H(s) is then substituted in Eq. (|5l|). In practice, the calculation can 
be simplified by using the perturbation theory value for O(s) in Eq. (|50|). We 
have argued above that perturbation theory is not valid for O(s). Nevertheless, 
using the perturbation theory value for Q(s) in Eq. (|50|) is justified because n(s) 
in perturbation theory has the same analytic structure in QCD. Thus using the 



perturbation theory value of II(s) in Eq. (51) is equivalent to using the perturba- 
tion theory value for Cl(s) in Eq. (|50|), even though the perturbative computation 
of O(s) is not valid. 

The OPE for II(s) is closely related to that for e + e~ — * hadrons, which depends 
on the time-ordered product of two electromagnetic currents. 

JL( S ) = Ci(p,s,a a (p))(Oi(tJL)), (52) 

where Cj are the coefficient functions, and Oi are the local operators. Since 
the contour C2 is a circle of radius m 2 in the complex s plane, one expects 
that logarithms in the uncalculated higher order corrections are minimized if one 
chooses \x = m T . The leading order operator is the unit operator. In the limit 
that the light quark masses are neglected, II^ vanishes, and we only need to 
compute Ht, giving 



2ttQ t (s) = \C T \ 2 {\V ud \ 2 + \V t 



2 (\ir |2 , 1 x y |2 

us 

1+ *s(Vs) +F3 
IT 




.(53) 



The first coefficient is the well-known result that the ratio <r(e + e — ► hadrons) /a(e + e 
qq) is 3(1 + a s (y/s)/ir). The next two coefficients are 



F 3 = 1.9857 - 0.1153n/, F 4 = -6.6368 - 1.2001n/ - 0.0052n). (54) 

Using the /3-function to write Q s (y / s) in terms of a s (m T ), evaluating the s integral, 
and setting nj = 3 gives |100| 1 



(55) 



R T = 3|C T | 2 (|K/ + |KJ 2 



1 + sV T > + 5.2023 v ' + 26.366 sK 1 ) + 



Assuming a s (m T ) rj 0.35, the series in brackets is 1 + 0.111 + 0.065 + 0.036 + . . ., 
so the terms are still numerically decreasing till order a 3 . One can use the value 
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of the last term as an estimate of the theoretical uncertainty in the perturbative 
value for the coefficient of the unit operator. This also tells us that we can neglect 
corrections from higher dimension operator that are smaller than about 3%. 

The 1/m^ corrections to R T arise from the quark mass corrections to the co- 
efficient of the unit operator in the OPE. At order m 2 /m^., the currents are no 
longer conserved, so one needs to compute both and Tlx- The only light quark 



mass contribution of any significance is the s-quark mass correction [100 1 , 



5R T = 3\V U . 



Hr(i + 



16a s (m T ) 
3^ 



(56) 



Using m s ~ 150 MeV as an estimate for the s-quark mass gives 5R T ~ —0.008, 
which is smaller than the error in the perturbation series. 

The and corrections in the OPE arise from the dimension four opera- 

tors F^ 1/ F >1U and mtpip, and the \jm\ corrections from the four-quark operators 
ipTipipTip, where T is some combination of 7 matrices. An analysis of these cor- 
rections, based on model estimates of the operator matrix elements indicates that 



these corrections are smaller than the uncertainty in the perturbation series [100]. 
The size of non-perturbative corrections can be determined directly from the ex- 
perimental data. Instead of considering the integral Eq. (^) that gives the total 
hadronic width, one compares the integral of dT/ds weighted with 

(1 — s / m%) k (s /m^) 1 with the corresponding moment of the experimental data. 
By studying the moments for different values of k and /, one finds that the non- 
perturbative corrections are about 3% [103, 104], and so are comparable in size 
to the uncertainty in the perturbation series. 

the 



The experimentally measured quantity is = 3.484 ±0.024 

ratio for r to decay into non-strange hadrons to the leptonic decay rate. This is 
given by Eq. (|55|), dropping V us , and gives 
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a s (m T ) = 0.34 ±0.03, 



(57) 



where we have assumed a theoretical uncertainty of 100% in the a 3 term. This 
value corresponds to 

a s {m z ) = 0.119 ±0.003, (58) 



7 a s FROM LATTICE GAUGE THEORY COMPUTATIONS 

The strong coupling constant a s can be determined from lattice gauge theory 
calculations of the hadronic spectrum. The basic procedure used is to choose a 
definition of a s , and measure its value on the lattice. One then has to set the scale 
at which a s takes on the measured value. The lattice scale can be normalized 
using the hadronic spectrum measured on the same lattice. Finally, one has to 
convert the lattice definition of a s to the value defined in the continuum in a 
scheme such as MS. 

There are several sources of systematic errors that limit the current accuracy 
in determining a s . Typically, a s is determined by determining the spectrum 
of heavy quark bound states on the lattice. There are corrections due to the 
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finite volume and finite lattice spacing a. The finite lattice spacing errors can be 
reduced by using improved actions, that are accurate to higher order in a. To 
some extent, one can estimate the error due to finite volume and finite lattice 
spacing by repeating the simulation on a larger lattice. The dominant systematic 
uncertainty is due to the quenched approximation, in which light quark loops 
are neglected. It is difficult to reliably estimate the systematic errors due to 
this approximation without doing a full simulation including dynamical fermions. 
Simulations with dynamical fermions are just starting to be done, and in a few 
years one should have more reliable estimates of a s . 

There is one important advantage to using a heavy quark system such as the 
T to determine a s . The leading correction to the T energy levels due to the light 
quark masses is linear in the quark masses, and can only depend on the flavor 
singlet combination m u + md + m s . Thus the light quark mass corrections can be 
computed to a good approximation using three light quarks of mass (m u + rrid + 
m s )/3. This avoids having to simulate almost massless dynamical quarks, which 
is very difficult. 

Lattice calculations can also be used to test theoretical calculations, and de- 
termine the regime in which perturbation theory is applicable. In the quenched 
approximation, one can study the scale dependence of the coupling constant on 
the lattice. This provides a check on the perturbation theory calculation with 
rif = 0. The result is in remarkable agreement with the perturbation theory 
result in the regime where the coupling constant is weak |105| , |106j ]. 

The Fermilab and SCRI groups use the S — P and IS* — 25 splittings in the T 
system to determine a s . There are some systematic deviations of the calculated 
numbers from their experimental values in the quenched approximation (nj = 0), 
which are dramatically reduced if one includes rif = 2 dynamical flavors. The 
value of a s (Mz) in the MS scheme is [107] 



a s (M z ) = 0.1159 ± 0.0019 ± 0.0013 ± 0.0019 (59) 

where the first error is due to discretization effects, relativistic corrections, and 
statistical errors, the second is due to dynamical fermions, and the third is from 
conversion uncertainties. 



More recent computations give (in the nj = 5 scheme) [ 108 ] 

a s {M z ) = 0.1174 ± 0.0024 (60) 

and POSH 



a s (M z ) = 0.1118 ± 0.0017. (61) 

An average of these newer values gives a s (Mz) = 0.115 ± 0.004, where we have 
included the difference between the two central values as an estimated additional 
systematic error. 



8 a s FROM HEAVY QUARK SYSTEMS 



Heavy quark bound states such as the T can also be used to extract a value for 
a s . If the bound state is treated using non-relativistic quantum mechanics, the 
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annihilation decays T — > fj, + fi~ , T — ► 555 and T — ► 755 can be computed as the 
product of the probability to find the quark-antiquark pair at the origin, times the 
annihilation rate for QQ at rest to decay to the final state. The relevant Feynman 
graphs are shown in Fig. |5[ The decay rate T — > ggg is the inclusive decay rate 
for T — > hadrons, and the decay rate for T — ► 7^ is that for T — > 7 + hadrons. 
The probability to find QQ at the origin, 1^(0)1 , is sensitive to the detailed 
dynamics of the QQ bound state. If one takes the ratio of decay rates, |^(0)| 
drops out, and the ratio of decay rates can be used to determine a s . 

The above qualitative discussion can be made precise using NRQCD (non- 
relativistic QCD) to calculate the properties of the T |110(| . NRQCD has an 
expansion in powers of v, the velocity of quarks in the bound state. The a s 
expansion is coupled to the v expansion, since v ~ a s in a Coulombic system. 
The NRQCD approach allows one to systematically factor the decay rate into 
short distance coefficients that are calculable in perturbation theory, and non- 
perturbative hadronic matrix elements that generalize the notion of the wavefunc- 
tion at the origin. The T wavefunction in NRQCD has different Fock components. 
The lowest order (in v) component is \QQ) and the first correction contains a 
gluon, \QQg). The \QQg) is referred to as the color-octet component, because 
the two quarks are in a color octet state. The NRQCD velocity counting rules 
show that the probability to find the T in the octet component is of order v 2 . 
The decay rate for T 



in NRQCD is [110| 



r(T -> M V 



2Im/ ee ( 3 5i^ 
M 2 



t o^Si) t) + 2Im ^ e 4 (35l \ r ViCS!) t 

where M& is the 6-quark pole mass, 0i( 3 £i) = ip'trx ■ X an d Vi{ 3 Si) = 
— i/A(ip^crx ■ x^cr(D) 2 -^ + h.c), and ip and \* annihilate quarks and antiquarks, 
respectively. The first term is the leading order contribution, and the second 
term is the v 2 correction. The coefficients f ee { 3 Si) and g ee ( 3 Si) can be computed 
from the first graph in Fig. || at lowest order in perturbation theory. The values 
are pO] 



2„2 



Im/ ee ( 3 5i 



to e 



1 



16^ 

7T 



Im.g ee ( 3 Si, 



■iira 2 e 2 
9 ' 



(63) 



where ej = —1/3 is the charge of the 6-quark, and the radiative correction to 
fee( 3 Si) has also been included. 



The decay rate for T — > hadrons is [110 ] 



r(T — > hadrons) 



2Im/ 1 ( 3 5 1 



Ml 



+r ( 8 )(T -» hadrons) 



T 



Oi( 3 Si) 



T> + 



2Im gl ( 3 5 1 ) 



T 



TM 3 ^) 



T 

(64) 



where is the contribution to the decay rate from the color-octet component 
of the T. In NRQCD, the color octet decay rate is v 4 suppressed relative to the 
color singlet decay rate. However, the color octet component can decay into two 
gluons, rather than three, so the color octet decay rate of order a 2 v 4 can compete 
with the relativistic correction to the color singlet decay rate of order a^v 2 . The 
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coefficient functions are [11C] 



Im/i( 3 5i 



Im 5 i( 3 5i) 



10 

243 

x 



(vr 2 - 9)a 3 s (M b ) 



a. 



1 + I -9.46(2)C> + 4.13(17)CU - 1.161(2)n/)-^ 
(19vr 2 - 132)a 3 (M fe ) 



1458 



(65) 



The a 3 term for Im/i( 3 Si) has the given value when the scale of the a 2 term is 
the 6-quark pole mass. 

The decay rate T(T — ► 7 + hadrons) has the form 



r(T — * 7+hadrons) 



2Im/ 7 ( 3 5i) 
Ml 



T 



Oi( 3 Si) 



T>+ 



2Img 7 ( 3 S 1 
Ml 



T 



V 1 ( 3 S 1 



T 



where the coefficient function is [ |110| 



(66) 



W 7 ( 3 Si) 
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(vr 2 -9)a 2 (M 6 )ae 2 



1 + -9.46(2)C F + 2.75(11)CU - 0.774(lW) — 



7T 



(67) 



The equations of motion can be used to relate the matrix elements of the 
S-wave and P-wave operators [111], 



T 



T 



M T - M b 
M b 



c?i( 3 5i) r) + o(v 



(68) 



The matrix element of the S-wave operator is the NRQCD analog of the wave- 
function at the origin in a potential model calculation. The matrix element is 
non-perturbative, and can be eliminated by considering ratios of decay rates. The 
matrix element of 7 7 i( 3 S'i) can be determined from Eq. ( |6"8| ) using estimates of 
the 6-quark pole mass M b , and the measured T mass, as was done in Ref. [111]. 
However, one can instead replace (My — 2M b )/M b by — 4/9a 2 , the lowest-order 
result for the binding energy for a Coulomb bound state. This reduces somewhat 
the uncertainty in the extraction of a s , since it eliminates any uncertainty from 
the pole mass. 

The experimental value of the ratio T(T — > hadrons)/T(Y — > £ + £~) = 39.11 ± 
0.4 ]112| , where £ = e,/j,,r gives a s (M b ) = 0.177 ± 0.01, using the ratio of the 
theoretical formula for the decay widths. The unknown octet decay rate has 
been estimated to be less than 9% in Ref. [Ill], and this has been included as 
a theoretical uncertainty. The decay rates Eq. (|62[)— ( 6€ ) have been written in 
terms of a s (M b ). One can instead rewrite them in terms of a s (/x) using Eq. (||), 
extract a s (fj,), and convert this into a s (M b ) using Eq. (g). We have include the 
uncertainty of a scale change by a factor of two in the theoretical estimate. The 
octet and scale uncertainties are comparable in size. 

The experimental value of the ratio T(T — ► 7 + hadrons)/T(Y — ► hadrons) = 
2.75 ±0.04 ±0.15 |113j ] can also be used to extract a s . It is convenient to use the 
experimental value of T(T — > hadrons)/T(Y — > £ + £~) to convert this to T(T — > 
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7 + hadrons)/r(Y — » = 1.075±0.06, before comparing with the theoretical 
results. This eliminates the theoretical uncertainty due to the octet component 
in the hadronic decay rate. The extracted value of a s (Mb) is 0.189 ± 0.01, where 
we have included a scale uncertainty as above. Averaging the two extractions 
gives a s (M b ) = 0.183 ± 0.01 which corresponds to a s {M z ) = 0.108 ± 0.004 



a s FROM HADRON-HADRON SCATTERING 



There are many process at high-energy hadron-hadron colliders which can con- 
strain the value of a s . All rely on the QCD improved parton model, and on the 
factorization theorems of QCD [114]. The rate for any process is expressed as 



a convolution of the partonic scattering amplitude of • and parton distribution 



functions discussed in section 0; see Eq 28 (note that here we use / rather than 



q as the sum on i,j runs over quarks and gluons. 

a = E/ dx^fi&M^fj&M^o^M) (69) 

The factorization scale M is arbitrary. As in the case of the scale \x used in a s (//) 
(see Eq [H] and the surrounding discussion), the exact result cannot depend on 
its choice. However as the processes er? ■ is only calculated to some finite order in 
perturbation theory, some residual M dependence will remain. As in the case of /U 
the sensitivity to M will be small if it is chosen to be a characteristic scale of the 
process; for example, in the case of the production of a pair of jets of momentum, 
pr, transverse to the direction defined by the incoming hadrons, M = px is a 
reasonable choice. 

The quantitative tests of QCD and the consequent extraction of a s which 
appears in of- are possible only if the process in question has been calculated 
beyond leading order in QCD perturbation theory. The production of hadrons 
with large transverse momentum in hadron-hadron collisions provides a direct 
probe of the scattering of quarks and gluons: qq — > qq, qg — > qg, gg — > gg, etc.. 
Here the leading order term in erf ^ is of order a? so the rates are sensitive to its 



value. Higher-order QCD calculations of the jet rates [115] and shapes are in 



impressive agreement with data |116|| . This agreement has led to the proposal 



that these data could be used to provide a determination of a s [117]. A set of 
structure functions is assumed and Tevatron collider data are fitted over a very 
large range of transverse momenta, to the QCD prediction for the underlying 
scattering process that depends on a s . The evolution of the coupling over this 
energy range (40 to 250 GeV) is therefore tested in the analysis. CDF obtains 
a s (M z ) = 0.1129 ± 0.0001 (stat.) ± 0.0085 (syst.) jn|]. Estimation of the 
theoretical errors is not straightforward. The structure functions used depend 
implicitly on a s and an iteration procedure must be used to obtain a consistent 
result; different sets of structure functions yield different correlations between the 
two values of a s . We estimate an uncertainty of ±0.005 from examining the fits. 
Ref. [117] estimates the error from unknown higher order QCD corrections to 



be ±0.005. Combining these then gives: a s (M z ) = 0.1129 ± 0.011 



27 



QCD corrections to Drell-Yan type cross sections (the production in hadron 
collisions by quark-antiquark annihilation of lepton pairs of invariant mass Q from 
virtual photons or of real W or Z bosons), are known [11£]. These processes are 
not very sensitive to a s as the leading piece in of ■ is of order a° . The production 
of W and Z bosons and photons at large transverse momentum begins at order 
a®. The leading-order QCD subprocesses are qq — ► 75 and qg — > jq. The next- 



to-leading-order QCD corrections are known |12C| ] [121] for photons, and for W/Z 
production [p.22||, and so an extraction of a s is possible in principle. 



Data exist on photon production from the CDF and D0 collaborations [123] 1 124 j 
and from fixed target experiments 125 ], Detailed comparisons with QCD pre- 
dictions |126| 1 may indicate an excess of the data over the theoretical prediction 



at low value of transverse momenta, although other authors [127] find smaller 
excesses. These differences indicate that while the process may be understood, 
no meaningful extraction of a s is possible. 

The UA2 collaboration ]128f | has extracted a value of a s (My/) = 0.123 ± 
0.018(stat.)±0.017(syst.) from the measured ratio R w = a(W + ljet)/a(W + Ojet). 
The result depends on the algorithm used to define a jet, and the dominant sys- 
tematic errors due to fragmentation and corrections for underlying events (the 
former causes jet energy to be lost, the latter causes it to be increased) are con- 
nected to the algorithm. The scale at which a s (M) is to be evaluated is not clear. 
A change from fj, = My/ to fj, = Mw /2 causes a shift of 0.01 in the extracted a s , 
and the quoted error should be increased to take this into account. There is also 
dependence on the parton distribution functions, and hence, a s appears explicitly 
in the formula for Rw, and implicitly in the distribution functions. Data from 
CDF and D0 on the W px distribution [12£] are in agreement with QCD but are 
not able to determine a s with sufficient precision to have any weight in a global 
average. 



The production rates of b quarks in pp have been used to determine a s IX 32|] 



The next to leading order QCD production processes [131] have been used. At 
order a s the production processes are gg — > bb and qq —* bb result in b-hadrons 
that are back to back in azimuth. By selecting events in this region the next-to 
leading order calculation can be used to compare rates to the measured value and 
a value of a s extracted. The errors are dominated by the measurement errors, 
the choice of \x and M, and uncertainties in the choice of structure functions. 
The last were estimated by varying the structure functions used. The result is 
a s (M z )= 0.1131ES 



10 CONCLUSION 

The previous sections have illustrated the large number of processes where quan- 
titative tests of QCD can be made and a value of a s extracted. Figure || shows 
the values of a s (Mz) deduced from the various processes shown above. The con- 
sistency and precision of these results is remarkable. Figure |7| shows the values 
of a s (p) and the values of [i where they are measured. This figure clearly shows 
the experimental evidence for the variation of a s (/i) with \i predicted by Eq.[|[ 
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An average of the values in Figure ^ and in Table [l] gives a s (M z ) = 0.1173, with 
a total x 2 °f 9 for twelve fitted points, showing good consistency among the data. 
The value from heavy quark systems contributes slightly more that one half of the 
total x 2 - If this result is omitted the average increases to 0.1185. All of the other 
results are within 1.1a of the average value. The error on the average, assuming 
that all of the errors in the contributing results are uncorrelated, is ±0.0014, and 
may be an underestimate. We have seen that in almost all of the cases discussed, 
the errors are dominated by systematic, usually theoretical errors. Only some 
of these, notably from the choice of scale, are correlated, it is important to 
note that the average is not dominated by a single measurement; there are many 
results with comparable small errors: from r decay, lattice gauge, theory deep 
inelastic scattering and the Z° width. We quote our average value as a s (Mz) = 
0.1173 ± 0.002, which corresponds to A^ 5 ) = 200±|| MeV using Eq. (jj. The 



reader may wish to consult other recent articles for different opinions |133| , 134]. 

Significant improvements in the precision in the near future are not likely. The 
accuracy of data from LEP will not improve. It is possible that a better under- 
standing of the jet rates in hadron-hadron colliders and a systematic treatment 
of the errors from the structure functions will lead to and improvement in the 
precision of the value of a s derived. In many cases where the data are quite 
precise, such as heavy quark system, theoretical uncertainties limit the precision. 
In the very long term precision at the 1% level may be achievable [135]. 
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Figure 1: Graphs that contribute to the propagator correction in the bubble chain 
approximation. 




Figure 2: Figure showing the kinematics of deep inelastic scattering ep — > eX 




Figure 5: Decay diagrams for T — > |U + /x , T — > and T — ► 7515. 
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Figure 7: Values of a s (fJ>) and the values of fj, where they are measured, these 
results are, in increasing value of fj,, r width, heavy quark decays, deep inelastic 
at \J < Q 2 > = 7 GeV, e + e~ annihilation rate at 35 GeV, e + e~ event shapes at 
58 GeV, the hadronic Z width, e + e~ event shapes at Z, 130 and at 189 GeV. 
The lines show our average and la errors. 



